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Neutrino losses from proto-neutron stars play a pivotal role to decide if these stars would be 
crushed into black holes or explode as supernovae. Recent observations of subluminous Type II-P 
supernovae (e.g., 2005cs, 2003gd, 1999br, 1997D) were able to rejuvenate the interest in 8-10 M0 
' stars which develop O+Ne+Mg cores. Simulation results of O+Ne+Mg cores show varying results 

in converting the collapse into an explosion. The neutrino energy loss rates are important input 
I parameters in core collapse simulations. Proton-neutron quasi-particle random phase approximation 

' (pn-QRPA) theory has been used for calculation of neutrino energy loss rates due to ^''Mg in stellar 

matter. The rates are presented on a detailed density-temperature grid suitable for simulation 
purposes. The calculated neutrino energy loss rates are enhanced up to more than one order of 
magnitude compared to the shell model calculations and favor a lower entropy for the core of these 
massive stars. 

PACS numbers: 23.40.Bw, 21.60.Jz, 26.30. Jk, 26.50.+X 

I. INTRODUCTION 

It is now more than forty years since Colgate & White [1] and Arnett [2] presented their classical work on energy 

• transport by neutrinos and antineutrinos in non-rotating massive stars. Despite immense technological advancements 
^ . . . . 

^ . in recent times, the explosion mechanism of core-collapse supernovae continues to challenge astrophysicists throughout 
the globe. When the iron core becomes unstable against the gravity force its inner portion undergoes homologous 
collapse whereas the outer portion collapses supersonically. The prompt shock that follows the bounce of the core 
stagnates (as it looses energy due to neutrino emission and endothermic dissociation of heavy nuclei falling through 
the shock )leading to a standing shock. Some additional source of energy is required to revitalize the postshock gas. 
One interesting mechanism to revive the shock was the "preheating" mechanism proposed by Haxton [3]. Haxton 
suggested that a significant fraction of the energy carried by the neutrinos could be utilized in preheating iron nuclei 
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outside the shock front rather than getting lost. This would lead to a reduction in the energy required from the shock 
wave to dissociate the infalling Fe, assisting the shock wave to retain more of its strength as it propagates through 
the iron core and thus increasing the chances of transforming the collapse into an explosion. However, the preheating 
mechanism continues to be debated and the part played by neutrinos in this scenario is far from being completely 
understood. Later, Bruenn and Haxton [4] worked on models simulating weak and strong shock cases and found 
out that in neither case is the energy transferred to the matter by neutrino-nucleus absorption significant in terms 
of preheating the infalling iron-like material. More recently, Langankc and collaborators [5] had some observations 
on the models of Bruenn and Haxton [4] and reported much larger preshock heating rates. They concluded that 
these rates act for too short a time to lead to consequences for shock propagation. Despite all this, the preheating 
mechanism in fact has been adopted by most coders, who are continually putting in more microphysics to make the 
simulations more realistic and there continues to be interest in this problem [6]. 

During the late stages of evolution, a star mainly looses energy through neutrinos and this process is fairly inde- 
pendent of the mass of the star. White dwarfs and supernovae, both have cooling rates largely dominated by neutrino 
production. Neutrinos are crucial to the life and afterlife of a supernova. It is natural to consider neutrino heating as 
a probable mechanism for shock revival as they dominate the energetics and dynamics of the post-bounce evolution. 
Having only weak interactions, they are nature's most efBcient means of cooling, transporting along around 99% of 
the released gravitational energy. As such it is no wonder that astrophysicists find it a formidable task to simulate a 
1% effect of explosion. 

White dwarfs located in a binary system may end their lives in two possible ways. They may accrete from a 
companion and achieve the Chandrasekhar mass thereby triggering a thermonuclear runaway of the object and 
ultimately exploding as a Type la supernova. In this case no remnant is left behind. Alternatively these massive 
white dwarfs, for sufficiently high mass accretion rates, may allow the formation of 0-|-Ne-|-Mg cores and due to the 
high prevailing central density (beyond 10^^ gcm~^) experience rapid electron capture that lead to the collapse of 
the core. This is termed as accretion-induced collapse (AIC). The end product is a neutron star (in some double- 
degenerate scenario, two white dwarfs in a short-period binary system may eventually coalesce to form a massive white 
dwarf that exceeds the Chandrasekhar mass limit and in this case transition to a black hole is possible if the total 
proto-neutron star mass exceeds the general-relativistic limit for gravitational stability [7]). The AIC of white dwarfs 
represents one instance where a neutrino mechanism leads undoubtedly to a successful, albeit weak, explosion [8] . The 
authors in Ref. [8] argued that the neutrino mechanism can successfully power explosions of low-mass progenitors and 
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AICs due to the limited mantle mass and steeply declining accretion rate. The ultimate fate of these white dwarfs 
are dependent on many factors, e.g. temperature of the environment, the mass accretion rate on the newly formed 
white dwarf, the mass of each partner white dwarf [7] and rapid electron capture rates [8] . The occurrence rate of the 
AIC of white dwarfs is not determined reliably, and are not expected to occur more than once per 20 - 50 standard 
Type la events (see for example [7]). 

Supernovae are also tipped in as probable sites for the r-process. The production site of the heavy r-process nuclei 
is associated with the accretion-induced collapse of an oxygen-ncon-magnesium (O+Ne+Mg) white dwarf in a binary 
system (e.g. [9]) or Type-II supernovae from 8-10 Mq (e.g. [10]). These 8-10 Mq stars form an electron-degenerate 
O+Ne+Mg core that does not undergo further nuclear burning. A few collapse simulators (e.g. Hillebrandt et al. [11]) 
demonstrated that the collapsing 0+Nc+Mg core explodes promptly. Others (e.g. Wooslcy and Baron [12]) showed 
that the shock generated at core bounce stalls rather than leading to a prompt explosion. Delayed explosion scenario 
of these cores was also reported (e.g. Mayle and Wilson [13]). Wheeler and collaborators [14] have suggested that the 
exploding O+Ne+Mg core could be a viable site for the r-process. However the question of whether O+Ne+Mg cores 
can explode hydrodynamically continues to be argued. For example, Gutierrez et al. [15] argued that the abundance 
of ^"^Mg was considerably reduced in updated evolutionary calculations. However the procedure adopted by them was 
not fully consistent as they kept the ratio of oxygen to neon constant while parameterizing the abundance of ^^Mg. 
Later Kitaura et al. [16] presented simulation result of these cores (keeping capture rates on ^^Mg as a key ingredient) 
using an improved neutrino transport treatment. Their outcome was not a prompt but a delayed explosion. 

A smaller iron core present at the onset of the core bounce as well as the smaller gravitational potential of the 
collapsing cores of a 8-10 Mq star do favor a prompt explosion. What remains certain is that microscopic and reliable 
electron capture rates and neutrino energy loss rates are needed for a careful study of the late stages of the stellar 
evolution of 8-10 Mq stars and can contribute effectively in the final outcome of the simulations of O+Ne+Mg cores 
on world's fastest supercomputers. 

Electron capture rates and the accompanying neutrino energy loss rates of matter in nuclear statistical equilibrium 
under conditions of high temperatures and densities are of prime importance in determining the equation of state of 
exploding stars. An accurate determination of neutrino emission rates is mandatory in order to perform a careful 
analysis of the final branches of star evolutionary tracks. A change in the cooling rates at the very last stages of 
massive star evolution could affect the evolutionary time scale and the configuration of iron core at the onset of the 
supernova explosion. 
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Compared to shell model calculations by Wildenthal et al. [17], the proton- neutron quasiparticle random phase 
approximation theory (pn-QRPA) gives similar accuracy in the description of light nuclei including beta-decay rates 
in sd-shell nuclide [18]. Nabi and Klapdor [19], later, calculated weak interaction rates for 709 nuclei with A = 18 to 
100 in stellar matter using the pn-QRPA theory. These included capture rates, decay rates, neutrino energy loss rates, 
probabilities of beta-delayed particle emissions and energy rate of these particle emissions (see also Refs. [20,21]). 
Since then these calculations were further refined with use of more efficient algorithms, incorporation of latest data 
from mass compilations and experimental values, and fine-tuning of model parameters (e.g. Refs. [22 25]). 

The evolution of the stars in the mass range 8-10 Mq develops central cores which are composed of ^^O, ^°Ne and 
^^Mg. Oda et al. [26] pointed out three different series of electron capture in the 0-|-Ne-|-Mg core of the 8-10 Mq 
stars and placed them in the order of low threshold energy as ^^Mg ^''Na ^*Ne, ^°Ne ^°0, and ^^O 

— >■ ^^N — > ^^C. The first series was regarded most important as the trigger nucleus, ^^Mg, has the lowest electron 
capture threshold. Oda and collaborators [26] actually renounced the last series in their calculations because of its 
high threshold energy which did not contribute significantly to the initiation of the collapse of the 0-|-Ne-|-Mg core of 
the 8-10 Mq stars. 

The Gamow- Teller strength distributions in ^^Mg and the calculation of electron capture rates on ^"'Mg using the 
pn-QRPA theory for 0-|-Ne-|-Mg core simulations were presented earlier by Nabi and Rahman [23]. Here I present 
for the first time the pn-QRPA calculations of the neutrino and antineutrino energy loss rates due to ^^Mg in stellar 
plasma on a detailed temperature-density grid suitable for the simulation codes of 0-|-Ne-|-Mg cores. In order to 
further increase the reliability of the calculated neutrino energy loss rates, I have incorporated all measured excitation 
energies (along with their log ft values) with definite spin and/or parity assignment, by either replacement (when 
they were within 0.5 MeV of each other) or manual insertion. The deformation parameter, S, is recently argued to 
be one of the most important parameters in QRPA calculations [27]. For the case of even-even nuclei, experimental 
deformations are available [28] and were employed in this work. All these steps were taken to ensure the most reliable 
calculation of pn-QRPA neutrino and antineutrino energy loss rates due to ^^Mg in presupernova and supernova 
environment. 

The paper is written as follows. Section II deals with the brief formalism of the pn-QRPA calculations. I present 
some of my calculated results in Section III. Comparisons with earlier calculations is also included in this section. I 
finally conclude in Section IV and at the end Table II presents the detailed calculation of neutrino and antineutrino 

energy loss rates due to ^''Mg suitable for simulation codes. 
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II. MODEL DESCRIPTION 

I used the pn-QRPA theory with a separable interaction to calculate the neutrino and antineutrino energy losses 
in stellar matter. The Hamiltonian of the problem was taken to be of the form: 

jjQRPA ^ jjsp ^ ypair + + yPP _ (1) 

For the single-particle Hamiltonian, H^p, the single particle energies as well as wave functions were calculated in the 
Nilsson model (which takes into account nuclear deformations). The BCS model was used to calculate the pairing force, 
ypoir rp^Q types of proton- neutron residual interactions were incorporated in these calculations namely the particle- 
hole, Vq^, and the particle-particle, Vq^, interaction. The two interactions (separable in form) were characterized 
by two interaction constants: x (for particle- hole interaction) and k (for particle-particle interaction). In this work, 
the values of x and k was taken as 0.001 MeV and 0.05 MeV, respectively. (The reader is referred to [18,29] for 
details of optimum selection of these parameters.) Other parameters required for the calculation of weak rates are the 
deformation, the pairing gaps, and the Q-value of the nuclear reactions. The calculated half-lives depend only weakly 
on the values of the pairing gaps [30]. Thus, the traditional choice of Ap = A„ = 12/ \/A{MeV) was applied in the 
present work. The deformation parameter is as an important parameter for QRPA calculations as pairing [27]. As 
such rather than using deformations from some theoretical mass models (as used in earlier calculations of pn-QRPA 
capture rates) the experimentally adopted value of the deformation parameters for ^^Mg, extracted by relating the 
measured energy of the first 2+ excited state with the quadrupole deformation, was taken from Raman et al. [28]. 
Q- values were taken from the recent mass compilation of Audi et al. [31]. 

An estimation of uncertainties of the model and theoretical errors is important to establish the reliability of the 
calculated rates. The uncertainties related to the pn-QRPA model was discussed in detail in Ref. [32]. There the 
accuracy of the pn-QRPA model was compared to the experimental data for both ^+ and /3~ directions (see Table 
I and Table II of Ref. [32]) and an estimation of theoretical errors was given. The decay rate of cvcn-cvcn nuclei 
generally exhibits more or less a smooth behavior with respect to the deformation parameter (see Fig.2 of Ref. [18] 
and discussions therein). The dependence of the half- life on the pairing gaps is small and rather smooth (see Fig. 
3 of Ref. [18]). The half-lives are not very sensitive to a change of the Gamow- Teller interaction strength x within 
a reasonable range and the relationship is almost linear [18]. On the other hand for large values of k the lowest 
eigenvalue becomes complex. This indicates the "collapse" of the QRPA model and it is crucial to avoid values of k 
which lie beyond the collapse. The redistribution of the /3 strength and the reduction of the sum of the Gamow- Teller 



6 

strength lead to a nonlinear dependence of the calculated half- lives as a function of k [29]. Figs. (7 - 10) of Ref. [29] 
show in detail how the observables of the pn-QRPA calculations vary with changes in k. These analysis do explain 
that a small and acceptable change in the QRPA parameters does not have a dramatic dependence on the calculated 
rates. 

The neutrino energy loss rates can occur through four different weak-interaction mediated channels: electron and 
positron emissions, and, continuum electron and positron captures. The neutrino energy loss rates were calculated 
using the formula 



\j - 



In 2 



(2) 



The value of D was taken to be 6295s [33] and the ratio of the axial vector to the vector coupling constant as -1.254 
[34]. B'^s arc the sum of reduced transition probabilities of the Fermi B(F) and Gamow- Teller transitions B(GT). 
The are the phase space integrals and are functions of stellar temperature (T), density (p) and Fermi energy 
of the electrons. They are explicitly given by 



u;2 _ \(w^ - wyF{±Z,w){\ - G^)dw, 



(3) 



and by 



= / w^v? - l(w„ wfF(±Z,w)G^dw. 

J Wl 



(4) 



In above equation w is the total energy of the electron including its rest mass, wi is the total capture threshold energy 
(rest-|-kinetic) for positron (or electron) capture. F(± Z,w) are the Fermi functions and were calculated according to 
the procedure adopted by Gove and Martin [35]. G± is the Fermi-Dirac distribution function for positrons (electrons). 

'E + 2 + Ef 



G+ = 



exp 



G- = 



exp 



kT 



kT 



1 



(5) 



(6) 



here E is the kinetic energy of the electrons and k is the Boltzmann constant. 

For the decay channel Eqt. (3) was used for the calculation of phase space integrals. Upper signs were used for 
the case of electron emissions and lower signs for the case of positron emissions. Regarding the capture channels, I 
used Eqt. (4) for the phase space integrals keeping upper signs for continuum electron captures and lower signs for 
continuum positron captures. 
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The total neutrino energy loss rate per unit time per nucleus is given by 

A^ = ^P,Aj;-, (7) 

ij 

where X"j is the sum of the electron capture and positron decay rates for the transition i ^ j and Pj is the probability 
of occupation of parent excited states which follows the normal Boltzmann distribution. 

On the other hand the total antineutrino energy loss rate per unit time per nucleus is given by 

Y = Y,PiKv (8) 

ij 

where X"^ is the sum of the positron capture and electron decay rates for the transition i ^ j. 

The pn-QRPA theory allows a microscopic state-by-state calculation of both sums present in Eqts. (7) and (8). In 
other words the pn-QRPA theory calculates the GT strength distribution of all excited states of parent nucleus in a 
microscopic fashion. This feature of the pn-QRPA model greatly increases the reliability of the calculated rates in 
stellar matter where there exists a finite probability of occupation of excited states. Further the pn-QRPA theory can 
handle any arbitrarily heavy system of nucleons since the calculation is performed in a luxurious model space of up 
to 7 major oscillator shells. 

Details of the calculations of phase space integrals and reduced transition probabilities can be found in [20]. 

III. RESULTS AND COMPARISON 

The Gamow- Teller strength distributions, B(GT±), for ^^Mg were calculated using the pn-QRPA theory. Quenching 
of the GT strength was taken into account and a standard quenching factor of 0.77 was used. The calculation was 
performed for a total of 132 excited states of ^^Mg up to excitation energies in the vicinity of 40 MeV (corresponding 
to first 100 excited states) in the daughters ^^Na and ^^Al. The total calculated B(GT+) and B(GT_) strength came 
out to be 3.34 (quenched) and Ikeda sum rule was satisfied in the calculations. 

In this section I present some of the results of the pn-QRPA calculated (anti)neutrino energy loss rates due to ^^Mg. 
The reported rates are also compared against previous calculations. 

Fig.l and Fig. 2 show four panels depicting the calculated neutrino and antineutrino energy loss rates, respectively, 
at selected temperature and density domain. It is pertinent to mention again that the neutrino energy loss rates 
(depicted in Fig.l) contain contributions due to electron capture and positron decay on ^^Mg whereas the antineutrino 
energy loss rates (Fig. 2) are calculated due to contributions from positron capture and electron decay on ^^Mg. The 
upper left panel, in both figures, shows the energy loss rates in low-density region {p[gcm^^] = 10°-^, 10^-^ and 10^-^), 
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the upper right in medium-low density region {p[gcm~^] = 10^'^, 10^'^ and 10^'^), the lower left in medium-high 
density region {p[gcm~^] = 10^-^, lO'"'^ and 10^'^) and finally the lower right panel depicts the calculated rates in high 
density region {p[gcm~^] — 10^'^, lO^*^'^ and 10^^). The (anti)neutrino energy loss rates are given in logarithmic scales 
(to base 10) in units of MeV.s~^. In the figures Tg gives the stellar temperature in units of 10^ K. One should note 
the order of magnitude differences in neutrino energy loss rates as the stellar temperature increases (Fig. 1). It can 
be seen from this figure that in the low density region the energy loss rates, as a function of stellar temperatures, are 
more or less superimposed on one another. This means that there is no appreciable change in the neutrino energy loss 
rates when increasing the density by an order of magnitude. There is a sharp exponential increase in the neutrino 
energy loss rates for the low, medium-low and medium-high density regions as the stellar temperature increases to Tg 
=5. Beyond this temperature the slope of the rates reduces drastically. One also observes that the neutrino energy 
loss rates are almost the same for the densities in the range (10 — lQ^)g/cm^ as a function of stellar temperature. For 
a given temperature the neutrino energy loss rates increase monotonically with increasing densities. 

As far as the antineutrino energy loss rates are concerned (Fig. 2) one again notes that the rates increase very 
sharply as Tg approaches 5. The rates are again almost superimposed on one another as a function of stellar densities. 
However as the stellar matter moves from the medium high density region to high density region these rates start 
to 'peel off' from one another. The antineutrino production rates are dominated by the capture of positrons due to 
^^Mg with a relatively smaller contribution coming from the beta-decay of ^^Mg (see Eqt. (8)). 

The current calculation was also compared with two earlier calculations of neutrino and antinciitrino energy loss 
rates due to ^''Mg. Fuller, Fowler and Newman [36] (hereafter FFN) compiled the experimental data and calculated 
neutrino and antineutrino energy loss rates (besides other weak-interaction mediated rates) for the nuclei in the mass 
range A = 21-60 for an extended grid of temperature and density. The GT strength and excitation energies were 
calculated using a zero-order shell model. For the discrete transitions, for which the ft values were not available, FFN 
took log ft = 5.0. Later Oda et al. [26] did an extensive calculation of stellar weak interaction rates of sd-shell nuclei 
in full (sd)"-shell model space. They used the effective interaction of Wildenthal [17] and the available experimental 
compilations for their calculations. 

Fig. 3 (Fig. 4) depicts the comparison of the pn-QRPA neutrino (antineutrino) energy loss rates with those calculated 
using shell model (OHMTS) [26] and those by FFN [36]. The ordinate represents the log (to base 10) of the calculated 
energy loss rates in units of MeV/s. The left panel depicts the different calculated energy loss rates at density 
lO'^gcm^^, the middle panel at density 10^gcm~^ and the right panel at density 10^^gcm~^. The shell model rates 
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are usually in good comparison with the corresponding FFN rates. The calculated pn-QRPA rates are enhanced. Table 
I shows the ratio of the calculated neutrino energy loss rates to those of Oda et al. [26] and FFN [36] calculations. In 
the table R^{QRPA/OHMTS) denotes the ratio of the reported neutrino energy loss rates to those calculated using 
the shell model whereas Ru{QRPA/FFN) gives the corresponding ratio for the FFN calculations. The reported rates 
are more than a factor 15 enhanced as compared to the shell model rates in the presupernova conditions. The shell 
model rates are slightly enhanced at Tg = 30 (at high densities the pn-QRPA neutrino energy loss rates again surpass 
the shell model rates). The collapse simulators should take note of this enhanced neutrino energy loss rates which favor 
cooler cores with a lower entropy. It is important to note from Table I that, in comparison to previous calculations, 
the QRPA rates are most enhanced around Tg ^ 1.5 — 2 (i.e. during the presupernova evolution of 0+Ne+Mg cores). 
The ratio then starts decreasing with increasing Tg. The enhancement of the QRPA rates around Tg ~ 1.5 — 2 may 
be traced back to the enhanced electron capture rates on ^^Mg [37] leading to an enhanced production of non-thermal 
neutrinos around these temperatures. It may also be seen from Table I that at high temperatures (Tg ~ 30) the ratios 
are in good agreement with previous calculations. 

The story is different for the comparison of antineutrino energy loss rates (Fig. 4). This time the shell model rates 
and FFN rates are much more enhanced compared to pn-QRPA calculations. Nevertheless the antineutrino energy 
loss rates are very small numbers and can change by orders of magnitude by a mere change of 0.5 MeV, or less, in 
parent or daughter excitation energies and are more indicative of the uncertainties present in the calculation of the 
excitation energies. 

Table II finally presents the calculated neutrino and antineutrino energy loss rates due to ^"^Mg on a detailed 
temperature-density grid suitable for simulation purposes. Here Column 1 shows the stellar density in logarithmic 
scales to base 10 (in units of gcm^^), Column 2 the stellar temperature in units of 10^ K. Stated also are the 
values of the Fermi energy of electrons in units of MeV in Column 3 whereas Column 4 and Column 5 display the 
corresponding calculated neutrino and antineutrino energy loss rates, respectively, in logarithmic scales to base 10 (in 
units of MeV.sec~^). The ASCII file of Table II is also available and can be received from the author upon request. 

IV. CONCLUSIONS 

There exists a wide variety of discrepant results in the supernova simulations of O+Ne+Mg cores including prompt 
explosions, delayed explosions and no explosions. Varying results of the simulations of 0-|-Ne-|-Mg cores ask for a more 
careful analysis of the precollapse evolution. A large number of input parameters are required for these simulations 
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to run and a reliable calculation of these parameters can certainly reduce the large uncertainties involved. Neutrino 
energy loss rates are one of the key input nuclear physics parameters responsible for cooling the stellar cores and 
thereby reducing its entropy. The pn-QRPA calculations of the neutrino energy loss rates due to ^'*Mg are presented 
here on a detailed temperature-density grid suitable for simulation codes. The pn-QRPA calculations gives similar 
accuracy in reproducing weak interaction mediated rates in sd-shell nuclide [18] and is employed here using maximum 
possible experimental incorporations to further increase the reliability of the calculated numbers. 

Recent simulations of 0+Ne+Mg cores [15,16] employ the weak interaction rates of shell model [38,26], respectively. 
The reported neutrino energy loss rates are up to a factor of 15 enhanced as compared to shell model results. These 
enhanced numbers favor cooler cores with lower entropies. Can the enhanced reported neutrino energy loss rates lead 
to any change in the outcome of simulation results on 0-|-Ne-|-Mg cores? Let us recall that in the simulations neutrino 
heating is a favored mechanism for revival of the stalled shock wave. The reported neutrino energy loss rates due to 
^^Mg and other nuclei of astrophysical importance (see e.g. Ref. [32]) employed preferably in a multi-dimensional 
model, with neutrino transport included consistently throughout the entire mass, accounting also for a complete 
treatment of multidimensional convection and burning phases might lead to some interesting outcome and answers. 
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Table I: Ratios of QRPA calculated neutrino energy loss rates (due to ^^Mg) to those of OHMTS [26] and FFN 
[36]. First column gives the grid point, containing the log (to base 10) of stellar density and temperature in units of 
gcm^"^ and 10^ K, respectively, at which these rates were calculated. 



(log(pye), Tg ) R^{QRPA/OHMTS) R^{QRPA/FFN) 



(1,1) 


1.24E+01 


1.13E+01 


(1,1.5) 


1.48E+01 


1.14E+01 


(1,2) 


1.48E+01 


9.84E+00 


(1,3) 


1.13E+01 


7.01E+00 


(1,5) 


5.51E+00 


3.91E+00 


(1,10) 


1.69E+00 


1.95E+00 


(1,30) 


7.14E-01 


1.61E+00 


(3,1) 


1.24E+01 


1.13E+01 


(3,1.5) 


1.48E+01 


1.14E+01 


(3,2) 


1.48E+01 


9.84E+00 


(3,3) 


1.13E+01 


7.01E+00 


(3,5) 


5.52E+00 


3.92E+00 


(3,10) 


1.69E+00 


1.95E+00 


(3,30) 


7.16E-01 


1.61E+00 


(7,1) 


6.71E+00 


5.98E+00 


(7,1.5) 


1.27E+01 


9.68E+00 


(7,2) 


1.41E+01 


9.35E+00 


(7,3) 


1.14E+01 


7.05E+00 


(7,5) 


5.85E+00 


4.08E+00 


(7,10) 


1.71E+00 


1.96E+00 


(7,30) 


7.16E-01 


1.62E+00 


(11,1) 


3.86E+00 


3.80E+00 


(11,1.5) 


3.87E+00 


3.79E+00 


(11,2) 


3.86E+00 


3.78E+00 


(11,3) 


3.84E+00 


3.67E+00 


(11,5) 


3.60E+00 


3.24E+00 


(11,10) 


2.90E+00 


2.20E+00 


(11,30) 


1.28E+00 


1.77E+00 
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Table II: Calculated neutrino and antincutrino energy loss rates due to ^ Mg for different selected densities and 
temperatures in stellar matter, log(pye) has units of gcm~^, where p is the baryon density and Yg is the ratio of the 
electron number to the baryon number. Temperatures (Tg) are measured in 10^ K. Ej is the total Fermi energy of 
electrons including the rest mass (MeV). are the total neutrino energy loss rates {MeVs~^) due to /3~^ decay and 
electron capture on ^''Mg . X" are the total antineutrino energy loss rates {MeVs~^) due to /3~ decay and positron 
capture on ^*Mg. The calculated rates are tabulated in logarithmic (to base 10) scale. In the table, -100 means that 
the rate is smaller than lO"^*^*^. 



logpYe 


^9 




A" 


A" 


logpYe 


^9 


Bf 


A" 


A" 




Tg 


Ef 


A" 




0.5 


0.5 


0.065 


-66.886 


-100 


1 


8-5 





-3-497 


-7.357 


2 


4.5 





-8-258 


-15.941 


0.5 


1 





-35.745 


-73.995 


1 


9 





-3.156 


-6.785 


2 


5 





-7.319 


-14.173 


0.5 


1.5 





-24.62 


-49.582 


1 


9.5 





-2.843 


-6.267 


2 


5.5 





-6.529 


-12.71 


0.5 


2 





-18.795 


-37.182 


1 


10 





-2.555 


-5.795 


2 


6 





-5-851 


-11-477 


0.5 


2.5 





-15.167 


-29.65 


1 


15 





-0.505 


-2.612 


2 


6.5 





-5.263 


-10.422 


0.5 


3 





-12.667 


-24.573 


1 


20 





0.773 


-0.8 


2 


7 





-4.745 


-9.506 


0.5 


3.5 





-10.825 


-20-907 


1 


25 


n 


1-697 


0-427 


2 


7.5 





-4-284 


-8-703 


0.5 


4 





-9.401 


-18.128 


1 


30 





2.416 


1.341 


2 


8 





-3.871 


-7.992 


0.5 


4.5 





-8.26 


-15.942 


1.5 


0.5 


0.162 


-65.908 


-100 


2 


8.5 





-3.496 


-7.356 


0.5 


5 





-7.321 


-14.175 


1-5 


1 


0-002 


-35-737 


-74.001 


2 


9 





-3-155 


-6-784 


0.5 


5.5 





-6.531 


-12.712 


1.5 


1.5 





-24.618 


-49.582 


2 


9.5 





-2.842 


-6.266 


0.5 


6 





-5.854 


-11.479 


1.5 


2 





-18.794 


-37.182 


2 


10 





-2.554 


-5.794 


0.5 


6.5 





-5.265 


-10-424 


1-5 


2-5 





-15-166 


-29-65 


2 


15 





-0-504 


-2-611 


0.5 


7 





-4.747 


-9.508 


1.5 


3 





-12.666 


-24.572 


2 


20 





0.774 


-0.799 


0.5 


7.5 





-4.287 


-8.705 


1.5 


3.5 





-10.824 


-20.906 


2 


25 





1.698 


0.429 


0.5 


8 





-3.873 


-7-994 


1-5 


4 





-9-399 


-18-126 


2 


30 





2-417 


1-343 


0.5 


8.5 





-3.499 


-7.359 


1.5 


4.5 





-8.259 


-15.941 


2.5 


0.5 


0.261 


-64.907 


-100 


0.5 


9 





-3.158 


-6.787 


1.5 


5 





-7.319 


-14.173 


2.5 


1 


0.015 


-35.669 


-74.07 


0.5 


9.5 





-2.845 


-6-269 


1.5 


5.5 





-6-529 


-12.711 


2.5 


1.5 


0-002 


-24-613 


-49-586 


0.5 


10 





-2.557 


-5.797 


1.5 


6 





-5.852 


-11.478 


2.5 


2 





-18.793 


-37.183 


0.5 


15 





-0.507 


-2.615 


1.5 


6.5 





-5.263 


-10.422 


2.5 


2.5 





-15.165 


-29.65 


0.5 


20 





0.77 


-0-803 


1-5 


7 





-4-745 


-9.506 


2.5 


3 





-12-666 


-24-572 


0.5 


25 





1.694 


0.424 


1.5 


7.5 





-4.284 


-8.703 


2.5 


3.5 





-10.823 


-20.906 


0.5 


30 





2.412 


1.337 


1.5 


8 





-3.871 


-7.992 


2.5 


4 





-9.399 


-18.126 


1 


0-5 


0-113 


-66.406 


-100 


1-5 


8-5 





-3-496 


-7.356 


2.5 


4-5 





-8.258 


-15-941 


1 


1 





-35.742 


-73.996 


1.5 


9 





-3.155 


-6.784 


2.5 


5 





-7.319 


-14.173 


1 


1.5 





-24.619 


-49.582 


1.5 


9.5 





-2.843 


-6.266 


2.5 


5.5 





-6.528 


-12.71 


1 


2 





-18.794 


-37-182 


1.5 


10 





-2.554 


-5.794 


2.5 


6 





-5-851 


-11.477 


1 


2.5 





-15.166 


-29.65 


1.5 


15 





-0.504 


-2.611 


2.5 


6.5 





-6.263 


-10.422 


1 


3 





-12.666 


-24.572 


1.5 


20 





0.774 


-0.799 


2.5 


7 





-4.745 


-9.506 


1 


3-5 





-10.824 


-20-906 


1-5 


25 





1-698 


0-428 


2.5 


7.5 





-4.284 


-8-703 


1 


4 





-9.4 


-18.127 


1.5 


30 





2.416 


1.342 


2.5 


8 





-3.87 


-7.991 


1 


4.5 





-8.259 


-15.941 


2 


0.5 


0.212 


-65.408 


-100 


2.5 


8.5 





-3.496 


-7.356 


1 


5 





-7.32 


-14.174 


2 


1 


0-005 


-35-72 


-74.018 


2.5 


9 





-3.155 


-6-784 


1 


5.5 





-6.529 


-12.711 


2 


1.5 





-24.617 


-49.583 


2.5 


9.5 





-2.842 


-6.266 


1 


6 





-5.852 


-11.478 


2 


2 





-18.793 


-37.182 


2.5 


10 





-2.554 


-5.794 


1 


6-5 





-5.264 


-10-422 


2 


2-5 





-15-166 


-29.65 


2.5 


15 





-0-504 


-2-611 


1 


7 





-4.746 


-9.507 


2 


3 





-12.666 


-24.572 


2.5 


20 





0.774 


-0.799 


1 


7.5 





-4.285 


-8.704 


2 


3.5 





-10.824 


-20.906 


2.5 


25 





1.698 


0.429 


1 


8 





-3.871 


-7.992 


2 


4 





-9.399 


-18.126 


2.5 


30 





2.417 


1.343 



15 







Ej 


A" 


\^ 




^9 


Ej 


\^ 




ZogpVg 


^9 


Ef 


A" 


x" 


3 


0.5 


0.311 


-64.406 


-100 


4 


0.5 


0.411 


-63.394 


-100 


5 


0.5 


0.522 


-62.276 


-100 


3 


1 


0.046 


-35.514 


-74,224 


4 


1 


0,209 


-34,692 


-76,046 


5 


1 


0,413 


-33,661 


-76,078 


3 


1.5 


0.005 


-24.602 


-49.597 


4 


1.5 


0.047 


-24,461 


-49.739 


5 


1.5 


0.265 


-23.727 


-50,473 


3 


2 


0.001 


-18.79 


-37.185 


4 


2 


0.014 


-18.759 


-37.216 


5 


2 


0.128 


-18.47 


-37.505 


3 


2.5 


0.001 


-15.165 


-29,661 


4 


2.6 


0.006 


-16,153 


-29,662 




2,6 


0,062 


-15,041 


-29,775 


3 


3 





-12.665 


-24.572 


4 


3 


0.004 


-12.66 


-24.578 


5 


3 


0.035 


-12.606 


-24.631 


3 


3.5 





-10.823 


-20.906 


4 


3.5 


0.002 


-10.82 


-20.909 


5 


3.5 


0.023 


-10.791 


-20.938 


3 


4 





-9.399 


-18,126 


4 


4 


0,002 


-9,397 


-18,128 




4 


0,016 


-9.379 


-18,146 


3 


4.5 





-8.258 


-15.941 


4 


4.5 


0.001 


-8.257 


-15.942 


5 


4.5 


0,012 


-8.245 


-15.954 


3 


5 





-7.319 


-14.173 


4 


5 


0.001 


-7.318 


-14.174 


5 


5 


0.009 


-7.31 


-14.182 


3 


5.6 





-6.528 


-12.71 


4 


5.6 


0.001 


-6,528 


-12.711 


6 


5.6 


0,007 


-6.622 


-12.717 


3 


6 





-5.851 


-11.477 


4 


6 


0.001 


-5.851 


-11.478 


5 


6 


0.006 


-5.846 


-11.482 


3 


6.5 





-5.263 


-10.422 


4 


6.5 


0.001 


-5.262 


-10.422 


5 


6.5 


0.005 


-5.259 


-10.425 


3 


7 





-4.745 


-9,506 


4 


7 





-4.744 


-9.506 


5 


7 


0,004 


-4.742 


-9,509 


3 


7.5 





-4.284 


-8.703 


4 


7.5 





-4.284 


-8.703 


5 


7.5 


0.004 


-4.281 


-8.705 


3 


8 





-3.87 


-7.991 


4 


8 





-3.87 


-7.991 


5 


8 


0.003 


-3.868 


-7.993 


3 


8-5 





-3.496 


-7.366 


4 


8,6 





-3,496 


-7.356 


5 


8,6 


0,003 


-3,494 


-7,357 


3 


9 





-3.165 


-6.784 


4 


9 





-3.166 


-6.784 


5 


9 


0.002 


-3.153 


-6,785 


3 


9.5 





-2.842 


-6.266 


4 


9.5 





-2.842 


-6.266 


5 


9.5 


0.002 


-2.841 


-6.267 


3 


10 





-2.664 


-6,794 


4 


10 





-2,664 


-5,794 


5 


10 


0.002 


-2.553 


-6,794 


3 


15 





-0.603 


-2.611 


4 


16 





-0.603 


-2.61 


5 


15 


0.001 


-0.503 


-2.611 


3 


20 





0.775 


-0.799 


4 


20 





0.775 


-0.799 


5 


20 





0.775 


-0.799 


3 


25 





1.699 


0,429 


4 


26 





1,699 


0,429 




25 





1,699 


0,429 


3 


30 





2.417 


1.343 


4 


30 





2.417 


1.343 


5 


30 





2.417 


1.343 


3.5 


0.5 


0.361 


-63.904 


-100 


4.5 


0.5 


0.464 


-62.865 


-100 


5.5 


0.5 


0.598 


-61.518 


-100 


3.5 


1 


0,115 


-35.167 


-74.571 


4.5 


1 


o,,sn9 


-34,187 


-76,552 


6.6 


1 


0,528 


-33.082 


-76.657 


3.5 


1.5 


0.015 


-24.668 


-49.632 


4.5 


1.5 


0,13 


-24,182 


-50.017 


5.5 


1.5 


0,423 


-23.198 


-51.003 


3.5 


2 


0.004 


-18.783 


-37.192 


4.5 


2 


0.044 


-18.683 


-37.292 


5.5 


2 


0.295 


-18.049 


-37.926 


3.5 


2.5 


0.002 


-15.162 


-29.653 


4.5 


2.5 


0.02 


-15.126 


-29.689 


5.5 


2.5 


0.18 


-14.803 


-30.013 


3.5 


3 


0.001 


-12.664 


-24.574 


4.5 


3 


0.011 


-12.647 


-24.591 


5.5 


3 


0.11 


-12.482 


-24.756 


3.5 


3.5 


0.001 


-10.822 


-20.907 


4.5 


3.5 


0.007 


-10.813 


-20.916 


5.5 


3.5 


0.072 


-10.72 


-21.009 


3.5 


4 


0.001 


-9.398 


-18.127 


4.5 


4 


0.005 


-9.393 


-18.132 


5.5 


4 


0.05 


-9.336 


-18.19 


3.5 


4.5 





-8.258 


-15.941 


4.5 


4.5 


0.004 


-8.254 


-15.945 


5.5 


4.5 


0.038 


-8.216 


-15.983 


3.5 


5 





-7.319 


-14.173 


4.5 


5 


0.003 


-7.316 


-14.176 


5.5 


5 


0.029 


-7.29 


-14.202 


3.5 


5.5 





-6.528 


-12.71 


4.5 


5.5 


0.002 


-6.526 


-12.712 


5.5 


5.5 


0.023 


-6.507 


-12.731 


3.5 


6 





-5.851 


-11.477 


4.5 


6 


0.002 


-5.85 


-11.479 


5.5 


6 


0.019 


-5.835 


-11.493 


3.5 


6.5 





-5.263 


-10.422 


4.5 


6.5 


0.002 


-5.261 


-10.423 


5.5 


6.5 


0.016 


-5.25 


-10.434 


3.5 


7 





-4.745 


-9.506 


4.5 


7 


0.001 


-4.744 


-9.507 


5.5 


7 


0.013 


-4.735 


-9.515 


3.5 


7.5 





-4.284 


-8.703 


4.5 


7.5 


0.001 


-4.283 


-8.703 


5.5 


7.5 


0.012 


-4.276 


-8.71 


3.5 


8 





-3.87 


-7.991 


4.5 


8 


0.001 


-3.87 


-7.992 


5.5 


8 


0.01 


-3.864 


-7.998 


3.5 


8.5 





-3.496 


-7.356 


4.5 


8.5 


0.001 


-3.495 


-7.356 


5.5 


8.5 


0.009 


-3.491 


-7.361 


3.5 


9 





-3.155 


-6.784 


4.5 


9 


0.001 


-3.154 


-6.784 


5.5 


9 


0.008 


-3.151 


-6.788 


3.5 


9.5 





-2.842 


-6.266 


4.5 


9.5 


0.001 


-2.842 


-6.266 


5.5 


9.5 


0.007 


-2.838 


-6.269 


3.5 


10 





-2.554 


-5.794 


4.5 


10 


0.001 


-2.553 


-5.794 


5.5 


10 


0.006 


-2.55 


-5.797 


3.5 


15 





-0.503 


-2.61 


4.5 


15 





-0.503 


-2.61 


5.5 


15 


0.003 


-0.502 


-2.611 


3.5 


20 





0.775 


-0.799 


4.5 


20 





0.775 


-0.799 


5.5 


20 


0.001 


0.775 


-0.799 


3.5 


25 





1.699 


0.429 


4.5 


25 





1.699 


0.429 


5.5 


25 


0.001 


1.699 


0.429 


3.5 


30 





2.417 


1.343 


4.5 


30 





2.417 


1.343 


5.5 


30 


0.001 


2.418 


1.343 



16 







Ej 


A" 


\^ 




^9 


Ej 


\^ 






^9 


Ef 


\^ 


x" 


6 


0.5 


0.713 


-60.357 


-100 


7 


0.5 


1.217 


-66.686 


-100 


8 


0.6 


2.444 


-46.023 


-100 


6 


1 


0.672 


-32-361 


-77.381 


7 


1 


1.2 


-29-976 


-80.038 


8 


1 


2-437 


-24-4 


-85-748 


6 


1.5 


0.604 


-22-592 


-51.61 


7 


1.5 


1-173 


-20.755 


-53.519 


8 


1.5 


2.424 


-16-885 


-57-677 


6 


2 


0.512 


-17.504 


-38.473 


7 


2 


1.133 


-15.966 


-40.037 


8 


2 


2.406 


-12.964 


-43.236 


6 


2-5 


0-405 


-14-351 


-30-465 


7 


2.5 


1.083 


-12-997 


-31-832 


8 


2.5 


2-383 


-10-511 


-34-45 


6 


3 


0.299 


-12-164 


-25-074 


7 


3 


1.021 


-10-959 


-26-287 


8 


3 


2-365 


-8-812 


-28.626 


6 


3.5 


0.214 


-10.515 


-21.214 


7 


3.5 


0.96 


-9.461 


-22.274 


8 


3.6 


2.322 


-7.663 


-24.248 


6 


4 


0.156 


-9.203 


-18-323 


7 


4 


0-871 


-8-307 


-19-223 


8 


4 


2-283 


-6.577 


-21-003 


6 


4.5 


0.117 


-8.127 


-16-072 


7 


4.5 


0-785 


-7-382 


-16-82 


8 


4.5 


2.24 


-5.792 


-18.449 


6 


5 


0.091 


-7.227 


-14.265 


7 


5 


0.698 


-6.618 


-14.877 


8 


5 


2.192 


-6.144 


-16.382 


6 


5.5 


0.073 


-6-461 


-12.777 


7 


5.5 


0.613 


-5-969 


-13-272 


8 


5-5 


2-139 


-4.597 


-14-67 


6 


6 


0.06 


-5-801 


-11.527 


7 


6 


0.534 


-5.404 


-11-926 


8 


6 


2-081 


-4.127 


-13.226 


6 


6.5 
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FIG. 1: (Color online) Neutrino energy loss rates due to ^*Mg, as a function of stellar temperatures, for different selected stellar 
densities . Densities are given in units of gcm~^. Temperatures are measured in 10^ K and logA^ represents the log (to base 
10) of neutrino energy loss rates in units of MeVsec~^ . 
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FIG. 2: (Color online) Same as Fig.l but for antineutrino energy loss rates. 
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FIG. 3: (Color online) Comparison of pn-QRPA calculated neutrino energy loss rates with those calculated using shell model 
(OHMTS) [26] and those calculated by FFN [36] as a function of stellar temperatures and densities. Ordinate represents the 
log (to base 10) of neutrino energy loss rates in units of MeVsec~^ . 




FIG. 4: (Color online) Same as Fig. 3 but for antineutrino energy loss rates. 




